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Abstract 

Data from HI and ZEUS on the structure and the quark and gluon densities 
of the proton are discussed. A brief excursion is made into the field of inclusive 
diffraction by deep inelastic scattering. The comparison of e~p and e^p scat- 
tering at large momentum transfers demonstrates clearly the presence of weak 
contributions in neutral current interactions. The comparison with the corre- 
sponding charged current results shows at Q"^ values above the masses squared 
of the heavy vector bosons the unification of electromagnetic and weak inter- 
actions. The new data are testing the validity of the Standard Model down to 
spatial resolutions of the order of 10~^^ cm. Intensive searches have been per- 
formed for a manifestation of instantons, for signs of compositeness of quarks 
and leptons and for the effect of extra dimensions. 



1 Introduction 

Ever since pointlike constituents, quarks, have been found in the nucleon the ques- 
tion has been raised whether quarks and leptons also have substructure. Obviously, 
the detection of quarks and leptons as extended objects would produce a revolution 

*This report is based on lectures presented at the IX Mexican School on Particles and Fields, 
Metepec, Puebla, Mexico, 9-19 August, 2000, and at the International School of Subnuclear Physics, 
Erice, Italy, 27 August - 5 September 2000. 
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in particle physics, e.g. by the possibihty of constructing quarks and leptons from 
common subconstituents. Deep inelastic lepton nucleon scattering is particularly well 
suited for the study of the spatial structure of quarks and leptons. The sensitivity to 
substructure depends on the virtuality of the exchanged current and, therefore, 
on the energy of the scattering partners. The highest energies are provided by HERA 
where by now a spatial resolution of 10~^^ cm or one thousandth of the proton radius 
has been achieved. The well defined initial state and the large energies and virtuali- 
ties make HERA an excellent place to probe the Standard Model also in other areas 
such as new currents and new particles. 

The first surprise encountered at HERA by HI p| and ZEUS in deep inelastic 
electron proton scattering has been a rapid rise of the proton structure function Fi 
as Bjorken-x tends to zero. The data are consistent with a power like rise, F2 oc a;~^, a 
behaviour which for a number of reasons , such as overlap of partons or a Froissart 
like bound for the total virtual photon proton cross section or the requirement of 
unitarity, should not persist down to infinitely small values of x. Instead, the parton 
densities should saturate at some x < xum- Saturation or not is a hot research topic at 
HERA. The discovery of large rapidity gap events in deep inelastic scattering [^, Q , 
which result from diffraction, has opened another avenue on this issue. The link 
provided by the optical theorem between the total cross section (which is proportional 
to F2) and the diffractive cross section suggests that the rise should be even faster in 
diffraction. 

This report presents new data from the two collider experiments HI and ZEUS at 
HERA with particular emphasis on the results for structure functions, diffraction and 
neutral and charged current processes at very large Q^. For a recent and comprehen- 
sive introduction to the physics at HERA see 0]. 



2 The HERA Collider and the Experiments 

The HERA collider can store electrons (positrons) of 30 GeV and protons of 920 GeV 
in two rings of 6.3 km circumference Table |l| lists some of the salient parameters 
of the machine 0. In order to maximize the luminosity up to 210 bunches of particles 
can be stored for each beam. The time interval between consecutive bunches is 96 
ns. The circulating e~ (e"*") beam becomes transversly polarized by the Sokolov- 



Ternov effect |]T0[- The measured specific luminosity is almost a factor of two larger 
compared to the design value. The maximum peak luminosity achieved so far is 
about 30% above the design value. The integrated luminosity per year provided by 
HERA for e-p collisions since 1992 is diplayed in Fig. ^ The total yearly luminosity 
increased by about a factor of two every year reaching 70 pb~^ in 2000. An upgrade 
program has been started which promises a factor of 3 - 5 increase in luminosity by 
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Table 1: HERA machine parameters. 



parameter 


electron 


ring 




proton ring 


circumference (m) 






6336 




beam energies, design (GeV) 




30 




820 


beam energies, achieved (GeV) 
e p cm. energy squared (GeV^) 




27.6 




920 






10^ 




magn. bending field (T) 




0.15 




5.25 


dipole bend, radius (m) 




610 




584 


max. circ. curr., design (mA) 




60 




160 


max. circ. curr., achiev.(mA) 




50 




110 


n. bunch buckets 




220 




220 


n. bunches (typical) 




180 




180 


time betw. cross, (ns) 






96 




max. lumi. design /achiev. 10^^cm~^s^^ 






1.5/2.0 




max. yearly lumi. deliv. per expt. (pb^^) 






70 




electron polarization 


50 - 


70% 







the insertion of additional magnets close to the interaction point. Operation in the 
new configuration will start in the year 2001. 

The data collected by HI and ZEUS between 1992 and 2000 for different beam con- 
ditions correspond to integrated luminosities of about 17 pb~^ for e~p and 116 pb~^ 
for e~^p collisions, per experiment. 



3 Structure Functions of the Proton 
3.1 Kinematics 

Inclusive deep inelastic scattering (DIS) by neutral current (NC) exchange, 

e{k) + p{P) e'{k') + anything 

can be described as a function of Bjorken-x and (see Fig. H). The basic quantities, 
in the absence of QED radiation, are: 

S = (1) 

Q' = -q' = -{e-er (2) 

X = (3) 
2P-q ^ ^ 

V = {q-P)/Mp (4) 
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HERA luminosity 1 992 - 2000 




Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 



Figure 1: Integrated yearly luminosity delivered by HERA per experiment during 
1992 - 2000. 

y - ^ (5) 

Q'^ ^ X - y- s (6) 
= ^!iiz^ + M2f^^fora;«l (7) 

X ^ X 

where e and e' are the four-momenta of the initial and final state electrons, P is the 
initial state proton four-momentum, Mp is the proton mass, s is the square of the 
ep cm. energy, —Q^ is the mass squared of the exchanged current, u is the energy 
transfer and y the fractional energy transfer from the incident electron to the proton 
as measured in the proton rest frame, y = i^max, ^max = s/ (2Mp), x is the fractional 
momentum of the proton carried by the struck quark and W is the 7*p cm. energy. 



4 




Figure 2: Diagram for deep inelastic ep scattering. 
3.2 Definition of the structure functions 

The differential cross section for deep inelastic scattering DIS can be expressed in 



terms of three structure functions, J-'2,^l,^3 Ijll 



[(1 + (1 - yY)r2 - y'J'L ± (1 - (1 - yr)xJ'3] (8) 



where a is the fine structure constant and the upper (lower) sign applies to e~ (e"*") p 
scattering. Since the contributions from JF^ and xJ^^ are expected to be small in the 
measured region, the radiatively corrected NC cross section can be written as: 



dxdQ^ xQ^ 

The structure function JF2 receives contributions from photon and Z° exchange and 
can be written as 

^2 = Fril + Sz). (10) 
where Fl'" denotes the contribution from photon exchange alone. 

The corrections 6z,l ,3 are functions of x and but are, to a good approximation, 
independent of JF2, i.e. they are insensitive to the parton density distributions. They 
were calculated from theory using structure functions which gave a good represen- 
tation of the data. In the measured region 6l is small except when y > 0.7 where 
6l ~ 0.12. The contributions from 6z,3 are negligible for < 1000 GeV^ and small 
up to ~ 5000 GeV^. Whenever, in the following, is below 1000 GeV^ the 
notations JF2, ^2"^™ and F2 are used interchangeably. 
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In QCD, ignoring Z° exchange, can be expressed in terms of the quark densities 
Q^) of the proton: 

^2 = Ee'a;g(x,g2), (11) 

where Cg is the electric charge of quark q and the summation is performed over all 
quarks and antiquarks. 

A recent in-depth discussion of the structure functions of the nucleon has been pre- 
sented in W\ . 



3.3 Structure function in the DIS regime 

The proton structure function T2 has been measured at HERA over a wide range in 
X and as shown in Fig. ^ At large x the HERA data overlap with those obtained 
by fixed target experiments. Until recently, most of the HERA results had come from 
data taken up to 1994, corresponding to about 2-3 pb~^ per experiment. HI and 
ZEUS have now presented preliminary analyses of data from 1996-7 which are based 
on an order of magnitude higher integrated luminosity. The region covered in x, 
has been enlarged substantially. At low ZEUS and HI have started to map out 
the transition region from photoproduction to deep inelastic scattering for x between 
10~^ and 10~^. At the high end, the structure function measurements have been 
extended up to ~ 3000O GeV^. Pro gress has also been made in the measurement 
of the charm contribution to F2 which provides for a direct test of the gluon density 
g[x^ Q^) extracted with QCD fits from the F2 data. 



The results on the x dependence of JF2 from HI and ZEUS |T^ are presented in 
Figs. ^, ^ for different intervals. The error bars show the statistical and systematic 
uncertainties added in quadrature. For < 100 GeV^ the typical statistical errors 
are 2% at low rising to 6% at ~ 100 GeV^. There is good agreement between 
the two HERA experiments. Also shown are the data from the fixed target experi- 
ments: BCDMS 0, E665 NMC [0 and SLAG [g which cover the region of 
^large^ x. In the region, where the HERA and the fixed target data overlap, good 
agreement is observed. 

Figure ^ shows the ^2 values as a function of for fixed x. For a; > 0.1 the data now 
span almost four decades in Q^. Scaling violations proportional to InQ^ are observed 
which decrease as x increases. 

The most striking feature of the HERA data is the rapid rise of JF2 as x — > which 
is seen to persist down to values as small as 1.5 GeV^. This rise accelerates with 
increasing as shown in Fig. |^ where ZEUS data from = 10, 22, 90 and 250 
GeV^ have been overlaid. Fitting ^2 data to the form J^2 = a + bx~^ leads to the 
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Figure 3: The x — plane: the regions covered by HI and ZEUS and by fixed 
target experiments. 



results given in Table [|T^ . The constant term a advocated in can be thought 
of representing the nonperturbative contributions such as expected from the aligned 
jet configuration, see [^|. The contribution from a decreases with Q^: for > 90 



GeV^ it is zero within the errors. The power A of the x dependent term is rather 
constant with while the coefficient b rises for between 10 and 90 GeV^ and 
appears to be driving the rapid rise of J^2- 

If the constant term is arising from a nonperturbative, hadron-like contribution then 
one may want to modify its contribution as follows: = ahad{^)°'''°"'^^^~^ + bx~'^, 

where ~ ^ is the square of the cm. energy of the virtual photon proton system, 
Wq = 140 GeV is the average W value of the data and apom(O) is the intercept 
of the pomeron trajectory obtained by fitting the cross section data from hadron- 
hadron scattering ||2^, ^ . The fit results for this form are presented in Table ^ The 



conclusions are the same as before. 



7 



ZEUS+Hl Preliminary 96/97 




Figure 4: Structure function from NC scattering as a function of x for fixed 
values of between 2.7 and 18 GeV^ as measured by HI, ZEUS. Also shown are the 
data from the fixed target experiments BCDMS, E665 and NMC The lines indicate 
QCD NLO fits to the data by HI and ZEUS. 
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ZEUS+Hl Preliminary 96/97 




Figure 5: Structure function from NC scattering as a function of x for fixed 
values of between 22 and 150 GeV^ as measured by HI, ZEUS. Also shown are the 
data from the fixed target experiments BCDMS, E665 and NMC The lines indicate 
QCD NLO fits to the data by HI and ZEUS. 
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• ZEUS 96/97 Preliminary 
■ HI 96/97 Prelim. ▲ HI 94/00 
A NMC, BCDMS, E665, SLAC 
ZEUS QCD Fit (Prelim.) 
HI QCD Fit (Prelim.) 

Ci(x) = 0.6(i(x)-0.4) 
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Figure 6: Structure function JF2 from NC scattering as a function of for fixed 
values of x as measured by HI, ZEUS. Also shown are the data from the fixed target 
experiments BCDMS, E665, NMC and SLAC. The dashed (sohd) lines indicate QCD 
NLO fits to the data by HI and ZEUS. 
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Figure 7: Structure function T2 from NC scattering as a function of x for fixed 
values of = 10, 22, 90 and 250 GeV^ as measured by ZEUS. For the curves see 
text. 



Fits performed without a constant term (a = 0) resuh in considerably larger values 
of /ndf for Q"^ < 40 GeV^. Thus, in this range of the data indicate the presence 



of a soft term a 7^ as has been noticed before 
rise with as shown in Fig. ^ see e.g. P^ , |25 



01. If a is set to zero, A is found to 



The curves shown in Figs. 
ZEUS 



I and § are the result of QCD NLO fits by HI [|T3 

m 



141 and 26, 27, 28, 29 based on DGLAP evolution 



The fits show that 

NLO DGLAP evolution can give a consistent description of the data over the full 
range. 
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Table 2: Parameters from the fit of ZEUS 1994 T2 data to the form JF2 = a + 6x ^ . 



parameter 


Q2=i0GeV2 


22GeV2 


90GeV2 


250GeV2 


a 


0.31±0.03 


0.14±0.05 


O.OliO.ll 


0.06±0.11 


b 


0.05±0.01 


0.10±0.02 


0.18±0.07 


0.11±0.06 


A 


0.36±0.02 


0.32±0.02 


0.32±0.05 


0.44±0.10 



Table 3: Parameters from the fit of ZEUS 1994 T2 data to the form T2 = 
ahad(w^)"''°'"^°^"^ + bx-^ with Wq = 140 GeV and apom(O) = 1.09. 



parameter 


Q2=iOGeV2 


22GeV2 


90GeV2 


250GeV2 


&had 

b 

A 


0.58±0.03 
0.013±0.003 
0.47±0.03 


0.32±0.06 
0.063±0.016 
0.37±0.03 


0.01±0.17 
0.18±0.08 
0.32±0.05 


0.11±0.15 
0.09±0.06 
0.47±0.11 



3.4 Structure function 

In the Quark Parton Model (QPM = zero order QCD), J-'l vanishes for spin 1/2 
partons. In LO QCD JF^ acquires a nonzero value due to the contribution from 
gluon radiation which is proportional to the strong coupling constant a^. A direct 
determination of J-'l requires the measurement of the DIS cross section at fixed x, 
for different values of y which can be accomplished e.g. by varying the ep cm. energy 
squared s. 

HI has shown that for a limited region of high y, Tl can be extracted from the T2 
measurements at a single value of s if these measurements are combined with a rather 



weak assumption on the validity of the DGLAP evolution [T^, |3T]]. At high y the 
factors 1 + (1 — y)^ and y^ which multiply T2 and jFj^, respectively, in the expression 
for the DIS cross section (see Eq. |^), are of comparable magnitude. With this in 
mind, the following procedure was chosen. The T2 values measured by HI for y < 
0.35 and by BCDMS at larger values of x are used to extract the parton distribution 
functions. The DGLAP equations allow to evolve the parton distribution functions in 
for fixed x and to predict JF2 at high y. Subtraction of the JF2 contribution yields 
then J-'l. Note, as shown above, NLO DGLAP gives a good description of the JF2 
data over four orders of magnitude in x and Q^, while for the determination of J-'l the 
evolution extends the maximum at fixed x by only a factor of two. Nevertheless, 
this analysis cannot strictly exclude the possibility that J^2 behaves differently than 
assumed. 

The longitudinal structure function J^l extracted by HI is shown in Fig. ^ as a 
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ZEUS 1995 
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0.5 
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ZEUSREGGE 
DL 



(GeV^) 



Figure 8: The power A as a function of obtained from fitting JF2 to the form 
jFo = bx~^. From ||25 . 



function of x. The full error bars represent the statistical and systematic uncertainties 
added in quadrature. The J-'l values are significantly above zero and a factor of 2 - 
3 below those of JF2. The dashed bands, which show J-'l as expected from the QCD 
NLO analysis, are consistent with the extracted J-'l values. 

This result lends also support to the procedure applied in the extraction of JF2 from 
the DIS cross ection where the (small) contribution from J-'l was taken from a QCD 
analysis, see above. 



3.5 Gluon density of the proton 

In the QPM (diagram(a) in Fig. |l^ the structure functions scale. Violations of 
scaling arise from QCD radiative effects which in LO (diagrams (b)-(d) in Fig. |T0]) 
have a simple mathematical expression: 
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Figure 9: Structure function JF^ from NC scattering as a function of x as determined 
by HI. The dashed band shows J^l as expected from the QCD NLO analysis. 



dJ^2 
rflnQ2 



Ee?^ J^%PAiy,Q')+PAix,Q% (12) 



Here, Pqq and Pqg are the quark and gluon sphtting functions and g{x, Q"^) is the gluon 
density of the proton. At small x, x < 10~^, the dominant contribution is expected to 



come from quark-pair creation by gluons (diagram (d) in Fig. |T0| and second term in 
Eq. |12D which offers the possibility to determine the density of gluons g in the proton 
rather directly. It is instructive to look first at the approximate relation derived by 



3^ who considers only the contribution from diagram(d): 



The scaling violations of F2, J^^q2 , measure directly the gluon momentum density 
xg{x,Q^). 
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q 7 q 

q \ \q 
(a) (b) 





Figure 10: DIS diagrams: (a) no QCD radiation, (b)-(d) lowest order QCD processes. 



HI and ZEUS determined the gluon density by performing a global DGLAP type 
QCD fit to the Ti data such as shown by the curves in Figs. ^ This method takes 
the contributions from the quark densities automatically into account but requires 
assumptions on the x dependence of the quark and gluon densities at the evolution 
scale Qo while for the approximate method these are not needed. 



The gluon momentum density xg{x,Q'^) of the proton as determined by HI at 
= 5, 20 and 200 GeV^ is shown in Fig. |ll|. Similar results have been reported by 
ZEUS [p5| , |26| . The precision now achieved for xg{x, Q"^) is around 8% at x = 5.10"^ 
and Q2 = 20 GeV^. 

It is instructive to convert the HERA data for JF2 and xg{x, Q'^) which give the quark 
plus antiquark and gluon momentum densities, respectively, into numbers of partons 
and compare these with those obtained at large x by |^ from the parton density 
set MRSDO'. The result is shown in Table ^ for = 20 GeV^ which corresponds 
to a spatial resolution of about 5.10"^^ cm. For high x, x > 0.06, one finds about 
2.4 quarks which is close to the canonical number of 3 quarks; in addition there 
are roughly 1.8 gluons. The parton numbers increase rapidly towards small x: for 
5.10""^ < X < 5.10-3 there are 4 times as many g, q and 15 times as many gluons as 
compared to the high- a; regime. 



Table 4: Equivalent number of partons of the proton for = 20 GeV^ determined 
at low X from the HI and ZEUS data and at large x from the MRSDO' set. 





X > 0.06 


5 ■ 10-^ < X < 5 ■ 10-=^ 




^ 2.4 
^ 1.8 


9 ± 1 
27 ± 2 
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Figure 11: The gluon momentum density of the proton as a function of x at Q"^ = 
5, 20, 200 GeV^ as determined by HI. 

3.6 Charm contribution to F2 

The structure function F2 measures the momentum densities of quarks in the proton 
summed over all quark flavours. Charm quarks are expected to contribute via fusion 
of the virtual photon with a gluon from the proton (boson-gluon fusion) (diagram 
(d) of Fig. |T0[). Therefore pair production of charm quarks in DIS offers another way 
to measure the gluon density of the proton. 

The charm contribution to DIS was determined by detecting D* production in a lim- 
ited 1]^* ,p^* region |Q ^ |5B|, |^. The measured cross section a{ep eD*X) 



was extrapolated to the full range in i]^* ,pY* with the help of a QCD model [^. Us- 
ing the branching ratio for c ^ D* measured at LEP the cross section a{ep eccX) 
and from this the charm contribution F2{x, Q^) were determined. In Figs. |l^, |13| F2 
is displayed as a function of x for fixed values, and for fixed x as a function of 
Q^. A comparison with the F2 data shows that for > Q GeV^ charm contributes 
about 20 - 30% of F2. This is in broad agreement with the ratio of 4/10 expected for 
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a democratic sea assuming massless quarks and neglecting the b quark contribution. 

The curves in Fig. 0, |13| show the prediction for F| using the gluon density as 
obtained from the QCD fit to F2. They are in good agreement with the data providing 
an important test for the QCD analyses and fits performed by HI and ZEUS on F2. 



in the NLO DGLAP scheme 

■ HI NLO QCD • HI Prelim, □ ZEUS 
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Figure 12: The charm contribution F2{x, Q^) to the proton structure function F2 
for fixed as a function of x as measured by HI and ZEUS. The curves show the 
predictions from the HI NLO QCD fit to the 1996-7 F2 data. 



In Fig. IT3 the dependence of F2 is compared for x values between 5.10 ^ and 0.02. 



The bands show the predictions of GRV94HO ||3^ for a range of charm quark masses 
between 1.2 and 1.6 GeV. The agreement is impressive. 



3.7 The rise of as a: ^ in the Ught of theory and models 



17 



Fi in the NLO DGLAP scheme 



2 

HlPrel. 



10' 



A ZEUS 

-i-i-i-l r 



HI NLO QCD 



fit to 



T 



x=0 ,00005 (s4 ) 



s=0.00013 (g4'") 



x=0,0002 (s4 ) 

0003 (x4°) 




10^ T tIO^ 

[GeVI 



1 
1 

1 



Figure 13: The charm contribution F2{x, Q^) to the proton structure function F2 
for fixed x as a function of as measured by HI and ZEUS. The curves show the 
predictions from the HI NLO QCD fit to the 1996-7 F2 data. 

3.7.1 The rise of .F2 as a; ^ and the W dependence of (7*°^ 

Neglecting contributions from exchange, the DIS cross section can be expressed in 
terms of the fiux of virtual photons times the total cross section for virtual - photon 
proton scattering, 0"*°*^ |^^, which is written in terms of the cross sections for the 
scattering of transverse and longitudinal photons, 

<4(a;, Q2) = aT{x, Q2) + aL{x, Q^). (14) 

The cross section defined in this manner can be interpreted in a way similar to the 
case of the interaction of real photons provided the lifetime of the virtual photon 
is large compared to the interaction time (in the proton rest frame), which means 
X -C l/(2Mp-i?p), where Rp is the proton radius, Rp = 4 GeV~^ This requirement 
is well satisfied if x <^ 0.1. The expression for JF2 in terms of ax and is 

U-,Q') = ^^^'r%. (15) 
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Figure 14: The charm contribution F2{x, Q^) to the proton structure function F2 
for fixed x as a function of as measured by HI and ZEUS. The bands show the 
GRV94HO predictions for masses of the charm quark between 1.2 and 1.6 GeV. 



At small X the expression can be rewritten in terms of the virtual-photon proton cm. 
energy W, ^ Q'^/x leading to 



<4^^^2(W^,Q^) (16) 
Equation |16| was used by ZEUS to determine from the 1993 J^2 data cr* ** . Fig- 



1 p 

ure |15| shows for the 1994 data Q'^a^°^p as a function of W from 20 to 260 GeV for 
fixed between 15 and 70 GeV^ |^|. The data cluster around a narrow band rising 



almost linearly with W . A fit of the form Q a °*p = a + b-W^ gave the value of e ~ 0.9 



(a fit with a=0 yielded a smaller value, ea=o ~ 0.5, though with considerable larger 
X^/ndf). If cr*°p is described in terms of a pomeron trajectory, 0"*°^ oc (H^2)'^jp(°)~^, 
the intercept ajp(O) = 1 + | 1.45. 

The observed rise of 0"*°^ is in marked contrast to the behaviour of the total cross 
section for antiproton-proton scattering and for real photon-proton scattering which 
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Figure 15: Q"^ function of W for fixed between 15 and 70 GeV^ 



as determined from the ZEUS data. The curves are proportional to the total cross 
sections ct*°* and a. 



IP 



-tot 

pp ■ 



are compatible with 0"*°^ ~ W , see curves in Fig. ^ Thus, the rise of as x 



or that of (T**^ as W 



00 signals the presence of a new phenomenon. In QCD, the 
rise is the result of the strong increase of the number of partons at small x which 
blacken the proton and reduce its transparency for virtual photons of high Q^. 



3.7.2 DGLAP evolution 



The DGLAP equations describe the evolution of the parton densities with Q^. In 
order to solve these equations one must provide the parton densities as a function of 
X at some reference scale Ql which should be large enough for perturbative QCD to 
be applicable Assuming a Regge-type behavior the small x dependence 

of the valence quark (v), sea quark (s) and gluon densities is of the form: 



xqy{x, Q"^) oc x"*^^"^ 
xqs{x, Q'^) oc x^~^>p 
xg{x, Q^) oc x^~'^F 



(17) 
(18) 
(19) 
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where and denote the intercepts of the reggeon and pomeron trajectories. For 
aji ~ 0.5 and ~ 1 one obtains xq^^x.Ql) oc x^'^ and xqs{x,Ql) oc xg{x,Ql) oc 
const. 

In the leading log approximation the x^'^ behavior of the valence distribution 
remains unchanged by the evolution while the sea quark and gluon distributions 
at small x become steeper. In fact, in perturbative QCD, JF2 is expected to grow 



faster than any power of ln(l/x) as x — [47, 48 1: 



X as{Q-'j 



exp^ 



144 In i ^^a^iQl) 



\ 33 - 2nf a,{Q^ 



(20) 



where ^/ is the number of quark flavors. The rise of JF2 as x — > can be accelerated by 
decreasing the reference scale Q^. This can be seen by applying Eq. ^ for a specific set 
of parameters, e.g. Uf = 3 and asiQ"^) = (^u^An^nQ'^/A'^ with A = 0.2 GeV. Starting 
with x-independent parton distributions at and parametrizing the J-'2{x, Q^) values 
obtained from Eq. |^ as J^2{x, Q^) = for 10"^ <x < 10'^ yields 



for Ql = A GeV2 at 
for Ql = 1 GeV2 at 



10 (20, 100)GeV2: A 
10 (20, 100)GeV2: A 



0.15 (0.21, 0.29) 
0.29 (0.32, 0.38). 



Hence, theory predicted that F2 would rise as x — but it could not predict the 
value of A and therefore the speed of the rise since the starting parton distributions 
at Qq were unknown. We will return to this point below in a discussion of the GRV 



model E9|, |50 



3.7.3 BFKL evolution 



The DGLAP scheme requires angular ordering and neglects terms proportional to 
In i, an approximation, which may run in difficulties as x ^ 0. The BFKL formal- 
|5T| does not impose angular ordering and resums terms proportional to In-. 



ism 



Based on the BFKL formalism which performs QCD evolution for fixed Q as func- 
tion of X the gluon density in the proton was predicted to rise as g{x, Q^) oc x 



-(l+A) 



as 



X — s> 0, where A ^ as(12/ ln2)/7r ^ 0.5 for Q"^ = 
NLO predict a considerably smaller value for A, 
fits to the data have been performed by . 



20 GeV^. BFKL-type calculations in 
A ^ 0.15 pi |S3|, g. BFKL inspired 
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3.7.4 Lopez- Yndurain model 



The authors of presented in 1980 - i.e. long before HERA came into operation 
- a NLO QCD model which predicted the rise of at small x observed at HERA 
with remarkable accuracy. According to the model, T2 should behave at small x 
as a power in x, T2 oc x'^'' where is independent of Q^, except for heavy flavor 
thresholds. Extending the scanty data available then down to x ~ 0.05 led to the 
prediction As = 0.37 ± 0.07. Adding a constant term to Ti a fit of the model to 
the new data from HI and ZEUS provided a good description of the measurements 
and yielded A^ = 0.355 ± 0.01 Dividing the data into different intervals 

indicated a possible but small rise of A^, with from 0.325 ± 0.01 at < 10 GeV^ 
to 0.355 ± 0.01 at > 100 GeV^. 

3.7.5 GRV model 

The rapid rise of JF2 at small x observed by the HERA experiments was anticipated 
in the GRV model In this model, the steep rise of JF2 at fixed and small 



X is generated dynamically by QCD evolution starting from an input scale Qo = fJ- 
through bremsstrahlung of gluons from quarks and by gluon annihilation into quark- 
antiquark pairs. At the evolution scale Qo = ji the input parton densities were chosen 
to be valence-like. 

The steepness of the rise depends on the input scale for which originally a value of 
/i^ = 0.34 GeV^ was chosen (GRV94). The GRV94 predictions gave a good account 
of the 1993 J^2 data from HERA but with the advent of the more precise T2 values 
from the 1994/5 and now from the 1996/7 data it became apparent that the rise 
predicted by GRV94 was somewhat steeper than measured. Recently, it was shown 



by the same authors ||50[ that agreement with the high precision data from HERA 
can be obtained by raising slightly the starting scale to /i^ = 0.4 GeV^ (GRV98). The 
predictions of GRV98 are shown in Fig.|16| together with 1994 and 1995 data from HI 
and ZEUS. For > 1.5 GeV^ the predictions follow closely the measurements. For 
lower the predictions fall below the data which is not surprising since there also 
nonperturbative contributions are expected. 

3.7.6 Haidt parametrisation 

Double-logarithmic scaling of J^2 with respect to x and has been investigated 
in ^]. The rise observed in the HERA data at small x is consistent with a 
logarithmic rise in x as well as in Q^. An economical parametrization which decribes 



the J^2 data of HI and ZEUS for x < 10"^ > q.h GeV^ has been obtained in §T 



J'2{x, Q') = m logio(l + %) logio(-) (21) 
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Figure 16: Structure function F2 from NC scattering as a function of x for fixed 
values of as measured by HI (1994/5 data) and ZEUS (1994 data). The F2 values 
are plotted as F2(x, Q"^) + i{Q'^) x 0.5 where for each set the value of i is indicated 
in the figure. The solid (dashed) lines show the GRV98 predictions calculated in NLO 
(LO). 

with Ql = 0.5 GeV^, xq = 0.04 and m = 0.41. 

Besides its simplicity, this parametrisation has the property that for fixed the 
total cross section for virtual photon proton scattering, atot{l*P) ~ ^-^2(2;, Q^), does 
not violate the Froissart bound for hadronic total cross sections. 

3.8 F2 in the transition region between photoproduction and 
DIS 

Virtual photon proton scattering at large with its rapidly rising cross section as 
W ^ 00 behaves markedly different from ordinary hadron hadron or real photon 
proton collisions. The rise is a sign for hard scattering as described by perturbative 
QCD (pQCD). This leads to the question where in does the transition from 
hadronic-type to hard scattering occur. When ^/Q^ is below the typical hadronic 
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transverse momenta of 200 - 500 MeV one expects confinement effects to dominate 
tfie interaction. On the otlier hand, the success of the GRV model may indicate that 
aheady at ~ 1 GeV^ hard scattering is the dominant type of interaction. In order 
to gain further insight, a precise mapping of the proton structure functions in the 
region from photoproduction up to of the order of a few GeV^ is essentiaL 

The structure function F2(x, Q^) vanishes as 0. This can be seen from its 

relation with the total virtual photon proton cross section, which is nonzero for 
= (Eq. [T^). Measurements of F2 at small and small x < 10"'' were published 
by HI m and ZEUS |2|, |3l. Figure shows F2 for fixed jC clS cL function of Q . 
A number of phenomenological models |6^, ^ |6^ have been put 



forward to describe the behaviour at low x and low Q^. The curves in Fig.|T^ show 



predictions of some of these models. The Vector Dominance type model DL p4| fails 



to reproduce the rise of JF2 for above ~ 0.5 GeV^. The BK model assumes 
a VDM-like component which dominates the region of low plus a hard QCD-like 
component for the high regime. The predictions are somwhat above the data at 
low Q^. The model CKMT ||6^ assumes at high the dominance of a bare pomeron 
with an intercept of ajp(O) ~ 1.24. At low the pomeron intercept is assumed to 
decrease due to rescattering corrections leading to 0^,(0) = 1.08, the value obtained 
from hadron-hadron scattering. The CKMT predictions are found to be below the 
data for < 0.6 — 1 GeV^. The GRV model considers only the hard scattering 
contribution. The GRV94 predictions for ^2 are close to zero for near the evolution 
scale Q2 = 0.34 GeV^. At = o.44 GeV^ GRV94 accounts for about 40% of the 
measured T2 and for about 80% at Q"^ = 0.57 GeV^. At = o.9 GeV^ basically 
all of the DIS cross section is attributed to hard scattering. The model ABY p7[ . 



which assumes a hard plus a soft component evolved in NLO-QCD, gives a rather 
good descripton of the full set of data. The comparison suggests that the transition 
from soft to hard scattering occurs at values somewhere between 0.8 and several 
GeV2. 



ZEUS studied also the scaling violations, dF2/d\nQ , in the transition region ||70 
25|, |63|. The logarithmic slope dF2/d\n.Q^ was derived from the data by fitting 



F2 = a + h\nQ'^ in bins of fixed x for W'^ ~ Q"^ /x > 10 GeV^. Figure H shows 
dF2/dliaQ'^ as a function of x (Caldwell plot). Also shown (on the top of the figure) 
for each x bin is the weighted mean of Q"^ {< >) which increases as x increases 
due to kinematics and detector acceptance. For x values down to 3 x 10~^ the slope 
dF2/dlnQ'^ increases as x increases. At lower x (equivalent to lower Q"^) the slope 
decreases. The prediction of the GRV94 model, for which dF2/dlnQ'^ was determined 
in the same manner as for the data, reproduces the data for a; > 3 x 10~^ (Q^ > 8 
GeV^). For smaller x the GRV94 slope keeps on rising while in the data it decreases. 

In order to gain further insight into the scaling violations at low x and ZEUS 
performed QCD NLO fits using all their data with 3 x 10~^ < x < 0.7 and > 1 



GeV together with those from NMC [71| and BCDMS |72]. A reasonable description 
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Figure 17: F2 as a function of for fixed values of x as measured by ZEUS, HI 
and E665. Also shown are predictions of various models, see text. 



of the data was achieved by the fits. Figure |T9| shows the singlet quark momentum 
density (xS = x Z]g[Q'(3^)+^(3;)]) and the momentum density of the gluon as a function 
of X for = 1, 7 and 20 GeV^. For >1 GeV^ the gluon density is much larger 
than the singlet quark density while at = 1 GeV^ the gluon density has become 
equal to or lies below the singlet quark density. Also, xS is seen to rise as x — > for 
all three values; xg , on the other hand, rises at = 7 and 20 GeV^ but may 
become constant at = 1 GeV^, or even zero. Such a behaviour was also found 

by m- 



One may tentatively conclude that at low Q^, of the order of 1 GeV^, the quark sea 
drives the gluon density while at higher the gluon drives the density of the sea 
quarks. However, this conclusion may be premature since the fits did not allow for 
a contribution from soft scattering which could still be substantial at > 1 GeV^ 
(see e.g. Fig. p!7| ). 
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Figure 18: dF2/dlnQ^ as a function of x calculated by fitting ZEUS F2 data in 
bins of X to the functional form a + blnQ"^. For each x bin the average value is 
indicated on top of the figure. The linked points labelled DL and GRV94 are from 
the Donnachie-Landshoff Regge fit and the GRV94 NLO QCD fit. In both cases, 
the points are obtained using the same range as for the experimental data. From 
ZEUS. 

4 Diffraction in deep inelastic scattering 



Diffraction has been studied extensively in hadron-hadron scattering at small mo- 
mentum transfers |7^. An elegant parametrization of the data has been provided 



by the Regge formalism through the introduction of a pomeron trajectory |75|. The 
hypothesis that diffraction may have a partonic component has been substanti- 
ated by the observation of high transverse energy jets produced in pp scattering |l77| . 
However, in hadron-hadron scattering both collision partners are extended objects 
which makes extraction of the underlying partonic process(es) difficult. In DIS, on 
the other hand, the virtual photon has a pointlike coupling to quarks. Hence, HERA 
offers a unique opportunity to study the partonic structure of diffraction since it 
gives, in a well defined manner, access to the regime of large photon virtualities and 
large energy transfers between the virtual photon and the target proton in its rest 
system, u = Q^/{2mpx) =2-20 TeV. 

Diffraction in virtual photon proton scattering has been studied at HERA in the 
quasielastic processes of vector meson production, 7*p Vp, where V = p^,u,(f), 
J/'^,T. While low mass V production (V = contributes more than 10% of 

the total cross section at = |^ it becomes negligible at large |j79|. However, 
diffractive dissociation of the virtual photon, 7*p XN (N= proton or a low mass 
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Figure 19: The quark singlet momentum distribution, xT,{x,Q'^) (shaded), and the 
gluon momentum density, xg{x,Q'^) (hatched), as a function of x at fixed values of 
Q2 = 1, 7, 20 GeV^ for two different evolution scales Q^^^. From ZEUS. 



nucleon system), into a large mass Mx, first recognized by the presence of a class of 
events with a large rapidity gap || remains a substantial fraction of the total DIS 
cross section also at large This has opened a window for a systematic study 

of diffraction in reactions initiated by a hard probe pO, pTl p2|, RSl Bl, [8^, 



4.1 t-dependence of the diffractive cross section 

The dependence of the diffractive cross section da-y*p^xp/ dMx on the square of the 
four-momentum transfer t between the incoming and outgoing proton was measured 
by ZEUS by detecting the scattered proton in the leading proton spectrometer (LPS) 



and the system X in the central detector |85]. The cross section is steeply falling 
with —t as shown in Fig. a fit of the form da^p^xp/dMx oc expipt) yielded h = 
7.2 ± 1.1 {stat)1lQ'Jj{syst) GeV~^. This shows that small momentum transfers between 
incoming and outgoing proton dominate, as expected for diffractive scattering. 
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Figure 20: The diffractive cross section da.y*p^xp/ d\A for events with a leading 
proton carrying more than 97% of the incoming proton momentum for 5 < < 20 
GeV2, 50 < < 270 GeV and 0.015 <5V(^i + <5^) < 0.5. From ZEUS. 



4.2 Diffractive structure function and cross section 



Isolation of diffractive events with the LPS is rather straightforward: detection of a 
proton scattered under very small angles and carrying a large fraction of the momen- 
tum of the incoming proton, xl = P^^''^ /Pp^^am ^ ^-^^ ensures a large rapidity gap 
between the outgoing proton and the system X. However, the event rate is limited 
by the acceptance of the LPS. 

In QCD, diffraction is characterized by the exchange of a colourless object, e.g. a 
colour singlet two-gluon system, between the incoming virtual photon and proton, 
producing a massive photon-like state and a nucleon, see Fig. |21|, or a low mass 
excited nucleonic system. In comparison with nondiffr active scattering, the exchange 
of a colourless system suppresses QCD radiation and therefore also the production of 
additional hadrons between the massive photon-like state and the nucleon or nucleonic 
system. Large acceptance for diffractive events has been achieved by requiring either a 
large rapidity gap between the nucleonic system N produced in the forward direction 
and the system X detected in the central detector, or by using the fact that in 
diffractive events most of the hadronic energy is carried away by the system which 
escapes detection leaving behind, in the region of the central detector, a low mass 
system. Therefore, by measuring the distribution of the mass of the hadronic system 
observed in the central detector the diffractive contribution can be separated from the 
nondiffractive one {Mx method). Analyses based on the first method were performed 



by HI 1 82, 84 1 and based on the Mx method by ZEUS E3, 86 
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Figure 21: QCD diagrams for diffraction in DIS. 



The results were presented in terms of the diffractive cross section 

da^*p^Xp/dMx and the diffractive structure function F:[^^^\x]p, f3, Q^) ||76|. The cross 



section for the process ep — > eXN is expressed in terms of the transverse (T) and 
longitudinal (L) cross sections, a^^-^ and crf^^, for ■y*p XN as: 

da'!;Ii^xN{Mx.W,Q') _ d{a^ff + aff) 



dMx dMx 

2tx Q2 daf"xN{Mx,W,Q^^ 



[22] 



a {l-yf + l dMxdlnW^dQ-^ 

The diffractive structure function of the proton can be related to the diffractive cross 
section in terms of the scaling variables xp ~ {Mx + Q'^)/{W^ + Q"^) and (3 ~ 
Q'^/{Mx + Q^)- In models where diffraction is described by the t-channel exchange 
of a system, for example the pomeron, x^p is the momentum fraction of the proton 
carried by the pomeron and (3 is the momentum fraction of the struck quark within 
the pomeron. One obtains |8l 



,a^UAM.,W,Q^) _ Z..,^ .rW..^,Q^ ,23) 



2Mx dMx (Q2 + iy2)2Q2 

For W"^ ^ Q^, Eq. ^ can be written as: 

1 dafIl^xN{Mx,W,Q^) An^a 
2Mx dMx g2(g2 + 



iP,x^,Q'). (24) 

If F^^^^ is interpreted in terms of quark densities then it specifies for a diffractive 
process the probability to find a quark carrying a momentum fraction x = jSXp of 
the proton momentum. 



4.3 Diffractive structure function measurement by HI 

HI presented their results for 7*p — > XN in terms of the diffractive structure function. 



The mass of the nucleon system N was restricted to Mn < 1.6 GeV. In Fig. 22 
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xipF^^^^ is shown as a function of xjp for fixed /? values and fixed between 4.5 
and 75 GeV^. The variation of xjpF^^^^ with (3 and is rather modest, indicating 
moderate scahng violations. In general, xjpF^^^^ is falling in the region xp < 10~^ 
followed sometimes by an increase at large xp values. 

The xjp dependence of ^2^^^'' is related to the W dependence of the diffractive cross 
section and, if analyzed in a Regge approach, to the Regge trajectories exchanged 
in the t-channel. By writing XjpF2^^\xp, (3,Q'^) = {C/Xp) ■ {xq/ XjpYF2^'^\l3,Q'^) 
one obtains n = 2(a^ — 1) if only the pomeron trajectory (here averaged over 
t) is contributing. Because of the rise of xjpF^^^^ seen at large xjp HI concluded 
that in addition to the pomeron a lower lying trajectory R is also contributing. The 
solid curves in Fig. 22 show the result of a two-component fit to the data. The 
dashed curves show the pomeron contribution alone as obtained from the fit. The 
fit yielded for the intercept of the pomeron trajectory a;^(0) = 1.203 ± 0.020(stat) ± 
{0.013{syst)~^Q^l^{model), a value which is above the results deduced from (soft) 
hadron-hadron scattering where ajp(O) = 1.08 and 1.096]'^o;oJ9 |2^ was found. 

HI fitted their data with a QCD motivated model, in which parton distributions 
are assigned to the leading and subleading exchanges. Figure |2^ shows the resulting 
contributions to the parton densities of the pomeron as a function of the fraction z 
of the pomeron momentum carried by the parton. Within this model the majority of 
the momentum of the pomeron is found to be carried by gluons. 



4.4 Diffractive cross section and structure function measured 
by ZEUS 

ZEUS determined the diffractive cross section and structure function for 7*p — > XN 



where M^r < 5.5 GeV |8^. The diffractive cross section is presented in Fig. |2^ as 
a function of W for various Mx and values. The diffractive cross section rises 
rapidly with W at all values for Mx up to 7.5 GeV. A fit to the form 

da:i^^xN{Mx,W,Q^) _ 

where a'^^^-f and the normalization constants h were treated as free parameters, yielded 
^diff ^ 0.507 ± 0.03A{stat)toMl{syst) which corresponds to a t-averaged IT^ = 1 + 
adiff/4 = 1.127 ± 0.009(stat)t[}:[J?^(s?/st). This value is consistent with the HI result 
since averaging over the t-distribution gives approximately cF^ = ajp(O) — 0.03. 

The diffractive cross section was compared with the measured total virtual-photon 
proton cross section. The ratio of the two cross sections, 

^ &dMxda'!^fl^xN/dMx 

^tot ' ^ ^ 
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Figure 22: The diffractive structure function xjpF<[^^^^ as a function of Xip for 
various P and values. The sohd curves show the results from the two-component 
Regge fit. The dashed curves show the pomeron contribution alone. From HI. 



is displayed in Fig. as a function of W for the different Mx bins and values. 



The data show that, for fixed Mx, contrary to naive expectations, the diffractive cross 
section possesses the same W dependence as the total cross section. The rapid rise of 
(Ttot with W, which is equivalent to the rapid rise of F2 as x ^ 0, in QCD is attributed 
to the evolution of partonic processes. The observation of similar W dependences for 
the total and diffractive cross sections suggests, therefore, that diffraction in DIS 
receives sizeable contributions from hard processes. The same W dependence for the 
diffractive and total cross sections was predicted in 



The diffractive structure function, multiplied by xjp is shown in Fig. ^ as a function 
of xip for different values of f3 and Q^. x ^F^ {x ^ , (3 , Q'^) decreases with increasing 
Xjp, which reflects the rapid increase of the diffractive cross section with rising W. The 
data are consistent with the assumption that the diffractive structure function F^^^^ 
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Figure 23: The sum of the hght quark and gluon distributions as a function of the 
momentum fraction z of the pomeron carried by the parton at different values of Q^. 
From HI. 

factorizes into a term depending only on and a structure function F^^"^^ which 
depends on {(3,Q^). The rise of x^F^^^^ with x^ can be described as XpF-^'^^^ oc 
with n = 0.253 ± .017 {stat)toml{syst). The data are also consistent with 
models which break factorization. 

Figure ^ shows F'^^'^\j3 , Q"^) = xoF^^^\xo, PjQ"^) where ^2^*-^'' was evaluated at 



xip = Xo = 0.0042. The data show that F-^^"^^ has a simple behaviour. For (3 < 0.6 
and < 14 GeV^, F^^'^^ is approximately independent of /3. For f3 < 0.8 also the 
data from different values are rather similar suggesting a leading twist behaviour 
characterized by a slow InQ^ type rescaling. For f3 > 0.9 the data show a decrease 
with (3 or Q"^. The approximate constancy of F^^"^^ for (3 < 0.9 combined with the 
rapid rise of ^2^''^^ as x^ decreases can be interpreted as evidence for a substantial 
partonic component in DIS diffraction dissociation. 

The behaviour of XjpF^^^\xjp, f], Q"^) is different from that of the proton structure 
function F2{x,Q'^), taken at x = x^, which rises gradually with Q^. It is in broad 
agreement with the conjecture that 



x^F^^^\x^,P,Q^) cx F2{x = x^,Q^)/\ogMVQl) where Ql = 0.55 GeV^. 
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Figure 24: 



as a 



The diffractive cross section da'^Ip^xN/^^x, Mm < 5.5 GeV, 
function of W for different Mx and values. The sohd curves show the result from 
fitting the diffractive cross section for each (VF, Q^) bin separately using the form 
da^f^^XN/dMx oc W-'-'' where a^^ff and the normalization constants were treated 
as free parameters. The dashed curves show the result from the fit where a'^*-^-^ 
assumed to be the same for all {W, Q"^) bins. From ZEUS. 



was 



4.5 Comparison with partonic models 



The data were compared with several partonic models of diffraction, NZ |8^, BPR |Q 
and BEKW |9l|| . Good agreement with the data can be achieved. The models provide 
a first glimpse of how the different components may build up the diffractive structure 
function. 

In the BEKW model basically three components build up the diffractive structure 



function, XjpF2 



D(3) 



Xjp, Q^) = ct ■ F^- + cl ■ + Cg ■ F^-^. The three terms repre- 
sent the contributions from transverse photons fluctuating into a gg or a qqg system 
and from longitudinal photons fluctuating into a qq system, and have the following 
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Figure 25: The ratio of the diffractive cross section, integrated over the Mx intervals 
indicated, a'^^^f — dMxCTylpi^xNy < 5.5 GeV, to the total cross section 

for virtual photon proton scattering, rf^l-^ = a^'^-^/cr*"^, as a function of W for the 
Mx intervals and values indicated. From ZEUS. 

functional dependence on x^p, /? and Q^: 

i^S=(-r^^^^ •/?(!-/?) (27) 

^.^,£^)„,«,.|.Hl + j|,)r/3=(l-2«= (28) 
= (^j'-f'' . + 0f (29) 

UtAQ') = 0-1 + 4l ■ ln[l + ln(^)] (30) 
In the model Xo,Ql, free parameters which have to be determined from the 
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Figure 26: The diffractive structure function of the proton multiphed by x^, 
function of Xp . The curves show the resuhs from the models of Nikolaev 
and Zhakarov (NZ), Bialas, Peschanski and Royon (BPR) and Bartels, Elhs, Kowalski 
and Wiisthoff (BEKW). 



data. The three terms F^, F^, F^-^ behave differently as a function of Q^. Except 
for a possible dependence of the power n^, does not depend on Q'^ . The term 
F^ is of higher twist but the power l/Q"^ is softened by a logarithmic factor; F^g 
grows logarithmically with Q^. 

It is instructive to compare the /3 and dependences of the three components. 
Figure p8| shows ctF^ (dashed), clF^ (dashed-dotted), CgF^g (dotted) and their sum 

XjpF^^^\xjp, f], Q"^) at Xjp = Xq (sohd curves) as a function of (3 for Q"^ = 8, 14, 27, 60 
GeV^. For /3 > 0.2 the colourless system couples predominantly to the quarks in 
the virtual photon. The region (3 > 0.8 is dominated by the contributions from 
longitudinal photons. The contribution from coupling of the colourless system to 
a qqg final state becomes important for (3 < 0.3. The last result is in contrast to 
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Figure 27: The structure function F^^'^^ (p, Q^) for YP XN, M^v < 5.5 GeV, for 
the values indicated, as a function of (3 as extracted from a fit to the measured 
XjpF^^^^ values. The curves show the fit results obtained with the BEKW model. 
From ZEUS. 

the HI observation (see above) that the large /3 region is dominated by the gluon 
contribution. Figure ^ shows the same quantities as a function of for (3 = 
0.1,0.5,0.9. The gluon term, which dominates at /? = 0.1 rises with while the 
quark term, which is important at /3 = 0.5 shows no evolution with Q^, i.e. is of 
leading twist. The contribution from longitudinal photons, which is higher twist and 
dominates at /5 = 0.9, decreases with Q^. 

In the BEKW model the -dependence of the quark and gluon contributions for 



transverse photons is expected to be dominated by the aligned jet configuration ^ 
and, therefore, to be close to that given by the soft pomeron. Writing oc 
(xq/xip)"''^ this implies rix ~ 2{ll^^°^^ — 1). However, perturbative admixtures in the 
diffractive final state are expected to have a somewhat stronger energy dependence, 
leading to an effective ut > 2(a^*°-^* — 1). The Xp dependence of the longitudinal 
contribution is driven by the square of the proton's gluon momentum density lead- 
ing to riL > riT- The fit of the BEKW model to the data indicates that transverse 
(longitudinal) photons dominate the region /? < 0.8 (/? > 0.8). 
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Figure 28: The three components {qq)T, ilig) and (gg)/, of the BEKW model 
building up the diffractive structure function of the proton and their sum F^^"^^ {fi^Q"^) 
as a function of (3 for = 8, 14, 27 and 60 GeV^, as obtained from a fit of the model 
to the data. From ZEUS. 



In a novel model (GBW) has been developed which links inclusive diffraction 
with the total cross section. In a frame where photon and proton are coUinear, the 
total 7*p and the diffractive cross sections can be written as [03, 03] : 



da^ii^ (t = 0) 1 



dt 



IGtt 



dh / d2:|^T,L(-2, r)no"(a;,r 



2\|2 



(31) 
(32) 



where '^T,L{z,'r) denotes the wave function for transverse (T) and longitudinal (L) 
photons, a{x, r^) the dipole cross section for the qq pair with the proton, z the 
momentum fraction of the photon carried by the quark and r the relative transverse 
separation between the quarks. 

The wave functions are determined by the photon-gg coupling and are known in 



QCD |9J]. The novelty of the model is an ansatz for the dipole cross section whose 
free parameters are determined by a comparison with the data for F2 ~ Q'^a^"*^. 
Given (T(x,r^), an absolute prediction can now be made for the diffractive cross 
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Figure 29: The three components {qq)T, ilig) and {qq)L of the BEKW model 
building up the diffractive structure function of the proton and their sum F^^'^\i3, Q^) 
as a function of for f3 = 0.1, 0.5 and 0.9 as obtained from a fit of the model to the 
data. The notation is the same as in the previous figure. From ZEUS. 



section. These predictions are found to give an almost quantitative representation of 
the ratio of diffractive cross section to total cross section as shown by Fig. [3^. This 
is a further step towards a quantitative description of DIS diffraction within QCD. 

An interesting discussion of DIS diffraction in terms of QCD radiation and the con- 



nection with the problem of confinement has recently been presented in p6 



5 Weak Interactions become strong: NC and CC 
Scattering at High 

In the standard model, electron proton scattering at low and medium (Q^ < 1000 
GeV^) proceeds almost exclusively through photon exchange. At higher values of 
substantial contributions are expected also from the exchange of the heavy vector 
bosons and Z°. The interference between photon and exchange contributes to 
e~p and e'^p scattering with opposite sign (see xF^ component in the expression for 
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Figure 30: Comparison of the predictions by the GBW model with the ZEUS data 
for the ratio of the diffractive to total cross section as a function of W, for the Mx 
intervals and values indicated. Figure taken from [ p5[] . 



the cross section below) which allows for a direct detection of the weak contribution 
in NC scattering. 

The structure functions can be expressed as sums over quark flavors of the proton's 
quark densities q{x, Q"^) weighted according to the gauge structure of the scattering 
amplitudes |Tl|, see also ||I2| . 



5.1 NC cross section 

For the neutral current (NC) reaction, e^p — > e^X, mediated by 7 and Z° exchange, 
the electroweak Born-level NC DIS differential cross-section can be written as 



^Vgo^.n(e^p) 



1^ 



(33) 
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where Y± = 1 ± (1 — y) . The structure functions and xF^ for longitudinally 
unpolarized beams may be described in leading order QCD as sums over the quark 
flavor f = u, ...,b of the product of electroweak quark couplings and quark momentum 
distributions in the proton 



xF. 



NC 



Y.xqj[Vl^A^-VfA 

f 



2 + ^vff + (4)2 + {Af) 



(34) 



where xqf = xqf{x,Q'^) ± xqf{x,Q'^) and xqf (xqj) are the quark (anti-quark) mo- 



mentum distributions. In leading order QCD, we have -F^^ 
and Af can be written as 



0. The functions Vi 



f 



L.R 
f 

aL,r 



Cf - {Ve ± ae) Vf 



(35) 



2ejsin26'w, are functions of the weak 



where the weak couplings, Oj = and Vi = Tf 
isospin, Tf = | (— |) for u, v {d, e), and the weak mixing angle, 6*^; Cj is the electric 
charge in units of the positron charge; and Xz is proportional to the ratio of Z°-boson 
and photon propagators 



4sin2^wCos2^w Q2 + M| ■ 



(36) 



The contribution of F^*^ to (fa^^^n/dxdQ'^ is predicted to be approximately 1.5% 
averaged over the kinematic range considered below. However, in the region of small 
X at the lower end of the range of the data shown below the F^ 
the cross-sections can be as large as 12%. 



contribution to 



5.2 CC cross section 



The electroweak Born cross section for the charged current reactions (see Fig. 31): 

e~p — >• fgX and e~^p — » (37) 

can be written as 



dx dQ"^ 



47rx VM^ + Q2 



nFCC(x, Q-)tY.xF^^{x, Q'')-y'Fi:^{x, Q^) 



i^CC 



2 l?CC/ 



where Gp is the Fermi constant and Mw is the mass of the W boson. 



(38) 
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Figure 31: Diagram for deep inelastic ep scattering by charged current exchange. 

The structure functions Fg^*" and xF^*^ , in leading-order (LO) QCD, measure sums 
and differences of quark and antiquark parton momentum distributions. For longitu- 
dinally unpolarized beams, e~p: 

F^"^ = x[u{x, Q^) + c{x, + Q') + s{x, Q')] (39) 
xF^^ = x[u{x, Q2) + c{x, Q2) - d{x, Q2) - s{x, Q^)] (40) 

F^"^ = x[d{x, Q2) + s{x, Q2) + u{x, Q2) + cix, Q2)] (41) 
xF^^ = x[d{x, Q^) + s{x, Q^) - u{x, Q^) - c(x, Q^)] (42) 

where u{x, Q^) is, for example, the number density of an up quark with momentum 
fraction x in the proton. Since the top quark mass is large and the off-diagonal 
elements of the CKM matrix are small, the contribution from the third generation 



quarks to the structure functions may be safely ignored [^. The chirality of the CC 
interaction is reflected by the factors Y± multiplying the structure functions. The 
longitudinal structure function, F^'^, is zero at leading order but is finite at next- 
to-leading-order (NLO) QCD. It gives a negligible contribution to the cross section 
except at y values close to 1, where it can be as large as 10%. 



5.3 Experimental results: NC scattering 



Both the HI and ZEUS experiments have previously reported cross section 
measurements, based on data collected in 1993, which established that the de- 
pendence of the CC cross section is consistent with the expectations from the W 
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propagator. The data from ZEUS demonstrated also that the CC and NC cross 
sections are of similar masr nitude for > M^. 



The high precision NC data presented recently by HI ||100| , |101|| (e p data from 1998- 
9 with 15 pb-^ e+p data from 1994-2000 with 46 pb'^) and ZEUS [|10g, |103l {cp 
data from 1998-9 with 16 pb~^; e^p data from 1996 with 30 pb^^) allow, for the first 
time, to see the contribution to the NC cross section. In Fig. ^ the e~p and 
e~^p cross sections are shown in terms of da/dQ"^ as a function of Q^. Both cross 
sections decrease by about six orders of mag nitude between = 500 and 40000 
GeV^, mainly governed by the photon propagator which leads to a behaviour of the 
form da/dQ"^ oc 1/Q^. For values above 3000 GeV^ there is a clear difference 
between the two charge states: the cross section for e~p scattering is larger than for 
e~^p, which demonstrates the presence of a weak contribution. 

The predictions of the Standard Model (solid curves) give a good description of the 
data. This is also true when the model is compared with the data for different values 



of X. This is shown in Figs. 33, 34 where the reduced cross sections 



_ XQ 1 daNC _ pNC I pNC /4o^ 

are given for fixed a; as a function of Q^. 

From the comparison of the e~p and e^p data the structure function xF^^ can be 
extracted. The dominant contribution comes from 7 interference, which is denoted 
by xF^^ . The reduced cross sections for the two charged states and the structure 
function xF^'^ as measured by HI | 100 | are presented in Figs. ^(a,b) for different 



intervals. 

In Fig. ^ xF^^ is shown as a function of x for values of 1500, 5000 and 12000 
GeV^. It is remarkable that only little dependence on is observed despite the large 
range. This is in agreement with QCD where a dependence on is expected 
only from scaling violations. Note also that xF^^ depends on the difference between 
quark and antiquark densities and is therefore primarily sensitive to the valence 
quark contribution. The data are compared with the results of a QCD fit performed 
previously by HI ||101| to their NC e~^p data taken in 1994-7 combined with the data 



from NMC and BCDMS (called HI 97 PDF fit). The fit certainly reproduces the 
qualitative features of the data. 

5.4 Experimental results: CC scattering 

The recent evaluations of CC cross sections by HI and ZEUS are based on an order 



of magnitude more data compared to the previous analyses, HI ||100|| (e p data from 
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Figure 32: The NC cross sections for e p and e^p as a function of as measured 
by HI and ZEUS. The curves show QCD-NLO predictions. 



1998-9 with 16 pb^^; e+p data from 1994-7 with 36 pb^^) and ZEUS 
data from 1998-9 with 16 pb"^ e+p data from 1994-7 with 48 pb"^). 

The combined CC cross section data from the two experiments are shown in Fig. |3^ 
for e^p and e~p scattering as a function of Q^. For < the cross sections 
show a slow decrease with which is mainly due to the shrinking phase space in 
X. The rapid fall at > is mainly driven by the propagator term +q2)2 - 
At low Q^, the cross sections for the two charge states are almost the same: the 
dominant contribution comes from scattering on sea quarks which is approximately 
flavor symmetric and contributes roughly equally in the two cases. At high the 
cross section for e~p scattering is larger than for e^p by almost an order of magnitude. 
Here valence quarks dominate the cross sections and different quark flavors contribute 
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Figure 33: The e~p NC reduced cross section for different values of x as a function 
of as measured by HI and ZEUS. The curves show QCD-NLO fits for 7 + Z'^ 
exchange. 



to the two charge states. Furthermore, the xF^ contribution is added to the e p cross 
section and subtracted for e^p. 



In Fig. ^ the reduced cross section, 



cr, 



cc 



2'KX 



+ Q' 



n 2 



<Pacc_ 
dxdQ'^ 



(44) 



as measured by ZEUS, is shown as a function of x for different intervals. As 
increases the reduced cross section for e~^p scattering becomes small with respect to 
the case of e~p scattering. 

While the charged vector bosons have been directly observed (CERN 1983) the 
effect of the W propagator on weak interaction cross sections, before HERA, had not 
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Figure 34: The e~^p NC reduced cross section for different values of x as a function 
of as measured by HI and ZEUS. The curves show QCD-NLO fits for 'j + 
exchange. 

been seen directly. As mentioned above, the fall-off of the charged current cross 
sections depends primarily on the W propagator. The new data from HERA allow a 
significant measurement of Mw thanks to the large range and the substantial size 
of the data samples. Using for Gp the measured value, the W mass was found to be: 



Mw = 
Mw 



79.9 ± 2.2{stat) ± 0.9{syst) ± 2.l{theor) GeV (HI) 
= SlAtll{stat) ± 2.Q{syst)tll{PDF) GeV (ZEUS). 



The values are in good agreement with the direct measurement of Mw 
0.056 GeV |Tn§. 



(45) 
80.419 ± 



45 



HI Neutral Current 



u r '\ 

, 1 _Q^=1500 GeV^ __Q^=5000 GeV" 



0.5 - 



• HI e p 

Vs=320 GeV - 



Hlep 

Vs=300 GeV - 




0.4 Q2=i500GeV^ ; Q^=5000 GeV^ 



0.2 - 




Q^=12000 GeV^_ 

— HI 97 
PDF Fit 



Q^=12000 GeV^ : 



Figure 35: The e~^p NC reduced cross section and xF^ for different values of as 
a function of x as measured by HI. The curves show the predictions from the HI 97 
PDF fit. 



5.5 Comparison of NC and CC cross sections 

In Figs. |3^, ^ the cross sections for NC and CC scattering in e~p and e~^p inter- 
actions are compared. At low Q^, where NC scattering arises predominantly from 
electromagnetic interactions (= photon exchange), the NC cross section exceeds by 
far the CC cross section. However, at values above 10000 GeV^ the two processes 
have about the same cross sections; at these large values of Q^, the weak force is of 
similar strength as the electromagnetic one: the HERA measurements show directly 
the unification of the electromagnetic and weak forces. 

It is instructive to compare this result with typical electromagnetic and weak particle 
decays. There, the weak force is about ten orders of magnitude smaller than the 
electromagnetic one. For instance, the decay time for the electromagnetic decay 
E° A7 is 7.4 ± 0.1 • 10"^° s while it is 4.1 ■ 10"^° s for the weak decay A ^ pvr", 
(note, both decays have similar cm. momentum, viz. 0.10 GeV and 0.074 GeV, 
respectively, and therefore similar phase space). 
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Figure 36: The structure function xF^^ for NC scattering as a function of x for 
different values of (from HI). The curves show the predictions from the HI 97 
PDF fit. 

5.6 Testing the Standard Model 

The large space covered in x and and their high precision make the NC and 
CC cross section measurements of HI and ZEUS a powerful testing ground of the 
Standard Model (SM). Uncertainties in the SM predictions arise from three sources: 
electroweak parameters, electroweak radiative corrections, and the parton momentum 
distributions including their higher order QCD corrections. A detailed discussion of 
these uncertainties can be found in ||10CI| , |101| , |102| , |104| 



For NC scattering, the electroweak parameters have been measured to high precision 
by other experiments and contribute less than 0.3% uncertainty for the SM pre- 
dictions. Higher order corrections for radiative corrections, vertex and propagator 
corrections and two-boson exchange are expected to be less than 1%. The primary 
source of uncertainties stem from the parton momentum distributions which have 
been determined from data mostly at low and then extrapolated to higher 
using DGLAP QCD evolution. The HERA data included in these fits are mostly 
from X < 10~^ and have little influence on the predictions for the high regime 
considered here. The parton densities (PDF) give a total uncertainty on the SM 
predictions for da^'^ /dQ"^ of 4% at = 400 GeV^ increasing to 8% at the highest 
Q"^ covered. 

For CC scattering, the main uncertainty of the SM predictions comes also from the 
uncertainties of the PDF's. The resulting uncertainty in da^'" /dQ"^ ranges from 4% 
at = 400 GeV^ to 10% at = loooo GeV^, and increases further at higher 
Q^. The large uncertainty at high is due to the ci-quark density which is poorly 
constrained at high x by the experimental data. 

In Fig. da^'" /dQ"^ as measured for e^p scattering is compared with the SM pre- 
dictions calculated with the CTEQ4D PDF's. Taking the errors of the data and the 
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Figure 37: The CC cross sections for e p and e+p as a function of as measured 
by HI and ZEUS. The curves show QCD NLO predictions. 
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Figure 38: The reduced CC cross sections for e p and e^p scattering as a function 
of X for different values (from ZEUS). The curves show the predictions using 
CTEQSD. 



uncertainties of the SM into account there is excellent agreement between the SM 
predictions and the data up to = 36000 GeV^. At higher Q"^ there are 2 events 
observed compared to 0.27 predicted. The corresponding data for CC scattering are 
presented in Fig. For < 5000 GeV^ the uncertainties in the PDF's - mainly 
due to the uncertainty of the d-quark contribution - prevent a precise test of the SM. 

The Standard Model (SM) provides a constraint between Gp and Mw- 

^ _ M| -[_ 

^ V2 {Ml - M^)M^ ^ ^ 

where Mz is the mass of the Z boson. The term contains the radiative corrections 
to the lowest order for Gp and is a function of aem and the masses of the fundamental 
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Figure 39: Comparison of the cross sections for e~p scattering by NC and CC 
exchange as a function of Q^, measured by HI and ZEUS. The curves show the 
predictions of the Standard ModeL 



bosons and fermions ||107|| . Using this constraint, assuming Mjj = 100 GeV and 
treating only Mw as a free parameter ZEUS found: 



Mw = 80. 50t°oitistat) to feisy St) ± 0. 31 {P D F) t o mi A Mt, AMh, AMz)GeV (47) 

The excellent agreement with the directly measured value of indicates that the 
SM gives a consistent description of a variety of phenomena over a wide range of 
energy scales. 



6 Search for new physics 
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Figure 40: Comparison of the cross sections for e^p scattering by NC and CC 
exchange as a function of Q^, measured by HI and ZEUS. The curves show the 
predictions of the Standard ModeL 



6.1 Inst ant ons 



The Standard Model contains a class of hard processes which cannot be treated 
perturbatively. These processes violate conservation of baryon and lepton number in 
electroweak interactions and of chirality in strong interactions ||108|| . Such processes 



are induced by instantons. QCD-instantons represent tunneling transitions between 
topologically different vacua [|109|| . Originally, instanton effects were expected to lead 
to observable effects only at energies of the order of 10^ GeV until it was realized 



that emission of gauge bosons can bring down the threshold ||110|, |111 



Up to now, experiments have been unable to demonstrate the presence of instanton 
effects. Deep inelastic lepton nucleon scattering (DIS) is particularly suited to search 
for such effects, firstly because DIS provides a hard scale (Q^), secondly, there exist 
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Figure 41: The e~^p NC DIS cross section for data and the Standard Model (SM) 
predictions (a) and the ratio of data to the SM prediction (b). The shaded region 
gives the uncertainty in the SM prediction due to the uncertainty in the PDF's. 



theoretical estimates of the instanton contribution involving quark-gluon 
boson) 



gauge 



interactions 
tial imj . 



112 1 and thirdly, the contribution is expected to be substan- 



First studies have been performed by HI ||114| , |115| , |116|| which placed limits on a 
possible instanton contribution. Preliminary results based on a much larger data 
sample (16 pb~^), taken in 1997, were reported recently by the same experiment ||117 |. 



The kinematics for DIS interactions involving instanton interactions is sketched in 
Fig.|43|. By instanton interaction a quark (momentum q', mass squared Q'"^) emitted 
by the virtual photon fuses with a gluon (momentum g = ^p) from the proton forming 
a hadronic system with cm. energy Wj. The cm. energy of the total hadronic system 
without the proton remnant is s. The instanton events have to be isolated from 
the standard DIS events depicted in Fig. H4l The main difference between the two 
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Figure 42: The e + p CC DIS cross section for data and the Standard Model (SM) 
predictions (a) and the ratio of data to the SM prediction (b). The shaded region 
gives the uncertainty in the SM prediction due to the uncertainty in the PDF's. 

classes of events lies in the fact that the particles produced by instanton interaction 
are phase space like distributed whereas those from standard DIS processes result 



from partons emitted in ladder-type diagrams, see Fig. 0. The challenge of the 
search for instanton effects lies in the task to devise suitable kinematic variables and 
find a kinematic region which maximize the instanton signal over the normal DIS 
background. 

The jet with the highest transverse energy {Etjet) is used to estimate the momentum 
q" of the current quark, see Fig.^. Knowing q" the momentum q' of the quark 
entering the instanton subprocess can be determined. The objects belonging to the 
jet are not used for computing the following variables. The iJ^-weighted rapidity 
mean fj = J2h Et^hVh/ J2h ^t,h of all remaining objects is used as an estimator of the 
instanton band which is defined by = ±1. The total scalar transverse energy Etb 
within this band is measured. Furthermore, the axis i is found for which the summed 
projections of the 4-momenta of all hadronic objects are minimal resp. maximal. 
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Figure 43: Diagram and kinematic variables for instanton interaction in DIS 



The relative difference between Ein = miriiJ^h \Ph ■ ^1 and Eout = maXiY^h \Ph ■ i\ is 
called Aft = {Bin — Eout) /Ein which measures the Et weighted azimuthal isotropy of 
an event, i.e. Af, = small (large) for isotropic (pencil-like) events. 

The search for instanton effects was performed with a sample of 275k DIS events. 
Figure 45 shows the distributions of Etjet, Eth and A^ for the data. They are well 
reproduced by the expectations for standard DIS as calculated with MEPS which 
incorporates the QCD matrix elements up to order as plus the emission of parton 
showers to all orders. A second model for standard DIS, CDM, where gluon emission 
is simulated via the colour dipole model, exceeds the data at large Etb and large Etjet- 
The expectations for the contributions from instanton processes as calculated with 
QCDINS are also shown ||118|] . Good discrimination between instanton and standard 
DIS contributions can be expected from Etf,: the instanton contribution is predicted 
to produce, on average, much higher Etb values. The predicted amount of instanton 
induced events is, however, at least a factor of about 20 below the expectations for 
standard DIS production determined with MEPS. In addition, at large Etf, there is 
a substantial uncertainty in the predictions for standard DIS contributions as shown 
by the difference between CDM and MEPS predictions. 

A summary of upper limits on instanton production from the HI analysis and the 
theoretical expectations for instanton production is presented in Fig. Instanton 
cross sections above 100 - 1000 nb are excluded. The instanton cross section predicted 
for x' > 0.35 and Q'^ > 113 GeV^ is below the upper limit derived from the data. In 
conclusion, an additional factor of 10 - 100 in discrimination against standard DIS is 
required in order to establish an instanton contribution of the predicted size. 
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Figure 44: DIS multi-ladder diagram. 



6.2 Deep inelastic scattering at high and high x 

During 1994-6 HI and ZEUS had collected sufficient data (14 pb'^ and 20 pb^\ 
respectively) for a ffist look at deep inelastic e^-p scattering beyond values of 10 
000 GeV^. Both experiments reported good agreement with the SM predictions for 
< 15 000 GeV^ and an excess of events in the region of > 15 000 GeV^, x > 
0.4 iTT9|, |l20i. 



HI observed 12 events with > 15 000 GeV^ where only 4.7 ± 0.76 were expected 
from SM calculations. Assuming the SM to be correct, the probability for such an 
excess was found to be about 1%. Further, these high events tended to cluster 
around x values of 0.4 - 0.5 corresponding to a positron-quark mass of M ^/xs ~ 
200 GeV which was suggestive for the production of a leptoquark or a R-parity 



violating SUSY state, see Fig. |^ In the region 187.5 < M < 212.5 GeV, y = 
Q^/{xs) > 0.4, 7 events were observed where 0.95 ± 0.18 were expected from SM 
processes. 

The ZEUS experiment found two events with > 35 000 GeV^ while 0.145 ± 0.013 
were expected from SM. For x > 0.55 and y > 0.25 four events were observed 
compared to 0.91 ±0.08 expected. In the SM the probability for such a fluctuation to 
occur in this kinematic region was 0.7%, and 8% for the entire region of > 5000 
GeV^ 

In 2000, both experiments presented preliminary analyses of the high Q^, high x 
region based on a factor 4-6 more e+p data and ffist high statistics data from e~p 
scattering. 
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Figure 45: Distributions of Etb, Af, and Etjet for the "instanton band" r/± 1. The 
data are compared with the expectations for standard DIS (MEPS, CDM) and for 
instanton production (QCDINS). From HI. 



The new resuhs from HI ||121|| are presented in Fig. ^ in terms of the electron- 



quark mass M, separately for the data taken in 1994-7 and 1999-0. While the earlier 
data show still some excess of events above the SM predictions (see histograms) near 
M = 200 GeV the data taken later agree well with the SM. Based on the data from 
1994-7 which correspond to an integrated luminosity of 37 pb~^ HI has placed tight 
constraints on the properties of first generation leptoquarks ||122|| . These studies were 
extended using the additional data taken during 1998-2000. The limits derived on the 
Yukawa coupling A for leptoquark (or R-parity violating squark) positron quark ver- 
tices based on the BRW model | 123|| are summarized in Fig. For masses between 



205 and 300 GeV these limits are more stringent than those from the TEVATRON 
or LEP. 



Figure. shows the distribution of events in the x-y plane as obtained by ZEUS ||124 
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Figure 46: Measured upper limit and expected cross section for instanton induced 
events as predicted by QCDINS as a function of x' and Q'"^. 



for e~^p data taken in 1994-7 and 1999-2000. Table |^ summarizes the number of events 
observed and expected at high Q^, and at high x, high y, including the results from 
e~p data taken in 1998-9. The combined data sets are in good agreement with the 
SM expectations and exclude a statistically significant excess. 

ZEUS has searched e^jet and V + jet final states produced in e^p interactions for 
leptoquarks and squarks [p.25| , |126|| . No signal has been observed. The limits derived 
for decays into e^jet or 17 + jet on the Yukawa couplings for scalar and vector LQ- 
lepton-quark and squark-lepton-quark couplings are summarized in Figs. H9|, pll. For 



vector states these limits are tighter than those obtained at the TEVATRON | |127|| and 
at LEP | 128|| . First generation scalar (vector) leptoquarks with an electromagnetic- 
type coupling strength, A = vivra = 0.31, are excluded for masses below 204 (265) 
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Figure 47: Diagram for leptoquark (or R-parity violating squark) production by ep 
scattering. 



Table 5: The observed and expected numbers of events with > 35 000 GeV^ and 
xda > 0.55,?/£)A > 0.25, respectively, from different years of HERA running, from 
ZEUS. 
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6.3 Search for lepton-flavor violation 



HI ||122|| and ZEUS [|129|| have searched for events of the type e^p — > IX where I is 
a final-state higher-generation lepton /i or r of high transverse momentum. Since no 
evidence was found for lepton-flavor violation the data were used to set constraints 
on couplings of leptoquarks (LQ) mediating lepton-fiavor-violating ep reactions. As 
an example. Table ^ shows the constraints derived for e ^ r transitions via fermion 
number F = LQ's. The data from HERA are shown together with the tightest 
limits from other experiments. For many transitions HI and ZEUS provide either the 
only or the smallest limit. 
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Figure 48: Mass spectra for NC DIS-like final states (points with error bars) and DIS 
SM predictions (histograms) for y > 0.1 and for a y cut which maximizes a possible 
signal from leptoquarks (squarks). Limits from HERA, LEP and TEVATRON. 

6.4 Substructure of quarks and electrons 

The SM assumes quarks and leptons to be pointlike. If, instead, they are extended 
objects the DIS cross sections at large will be reduced in comparison to the SM 
prediction. The modification resulting from formfactors of electron and quarks can 
be written as 



YFeiQ')F;m (45 



dQ^ dQ 



where '^'^^2 represents the SM cross section. Possible modifications of the qq — gluon 
vertex from extended quarks are expected to be small and are neglected. 



With the assumption that electrons are pointlike HI |130| has derived an upper 



limit for the quark radius of Rq < 1.7 ■ 10 cm. Assuming the electron and quark 
formfactors to be the same and of the form -Fe(Q^) = -Fg(Q^) = 1/(1 + |-R^Q^) 
the data presented by ZEUS [ |102|| for the ratio ^^/^Q2(g°j(,/j in e+p NC scattering 
lead Q to a limit of i?e,g < 1-5 ■ 10^^^ cm for the radii of electrons and quarks. The 
^Extracted by the author. 
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Figure 49: Exclusion limits from HERA, LEP and TEVATRON at 95% CL on the 
Yukawa coupling A as a function of the LQ mass for a scalar LQ with F = (top) 
and F = 2 (bottom) described by the BRW model. The grey and hatched domains 
are excluded. 
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Figure 50: The distribution in xda and i/da of e~^p data taken in 1994-7 (left) and 
1999-0 (right). The curves indicate constant values of Ql?^ = xdaUdas for Qf)^ = 
5000 and 35 000 GeV^. In the left distribution the open circles mark the five "high-x 
high-y" events from the 1994-6 data set. From ZEUS. 



HERA limits can be compared with the result Rq < 1 ■ 10^^^ cm obtained at the 
TEVATRON ||131|| from Drell-Yan pair production assuming pointlike leptons. An 
analysis of the contributions of anomalous magnetic dipole moments to the Zqq vertex 
yielded Rq < 1.2 ■ 10~^^ cm for u and d quarks 



6.5 Extra dimensions 



Starting from the Kaluza-Klein model [133] it has recently been suggested that the 
gravitational scale in 4 -|- n dimensions may be as low as the electroweak scale [|134| , 
|135|| . While in conventional string theory the extra dimensions are assumed to be 



compactified, in models with large extra dimensions the spin 2 graviton propagtaes 
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Figure 51: Exclusion limits at 95% CL on the Yukawa coupling A as a function of 
the LQ (squark) mass for scalar and vector states decaying into e^q and Vq channels, 
derived from jet and Vjet final states. The horizontal lines marked A = 0.31 
indicate an electromagnetic coupling strength A = \/47ra = 0.31. From ZEUS. 
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Figure 52: 95% CL rejection limits on "112 ^"^^ for F = LQs (in units of 10 ^ 

GeV~^). The first column indicates the quark generations coupling to LQ — e and 
LQ — T respectively. The second row indicates the different F = LQs. The ZEUS 
preliminary results are shown in the third line (bold) for each case. The process 
providing the most stringent limit on the given LQ is shown in the first line of each 
case together with its value (second line). If there is no constraint from low energy 
experiments, the ZEUS result is compared to the HI result. ZEUS limits which are 
comparable (within a factor of two) or better than the low energy results are enclosed 
in a box. The * indicates the cases where only a top quark would be involved. From 
ZEUS. 
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into the extra spatial dimensions. As a result the gravitational force maybe strong in 
4 + n dimensions and much reduced in our four dimensional world. In lepton-quark 
scattering this mechanism leads effectively to a contact interaction with a coupling 
coefficient of 7]g = A/M| pGj . 



HI ||137|| has searched for the presence of such a contact interaction in e p and e p 



NC interactions by studying the ratio '^^J^^2^sm\ ' -^^S- No evidence for contact 
interactions was observed. The 95% CL. lower limits are Ms = 0.63 TeV (A = +1) 
and Ms = 0.93 TeV (A = —1). Similar limits have been derived by LEP experiments 
from e~^e~ annihilation ||138 . 



6.6 FCNC and single top production 

In the SM, flavour- changing neutral current (FCNC) transitions are forbidden at 
the tree level, and a transition between a top quark and a charm or up quark can 
only occur via loops which renders the cross section extremely small. However, some 
models beyond the SM predict anomalous effective couplings of the top-quark of the 
type tuV and tcV with V = j, Z^, which could give rise to experimentally detectable 
signals [ |139|| . In ep collisions they can lead to single top production. Since the top 
decays almost exclusively into bW a promising event signature are topologies with 
an isolated lepton, missing transverse momentum and a large hadronic transverse 
momentum. 

Both HERA experiments have searched for such events. No candidate was observed 
by ZEUS in e^p data from a total of 82 pb~^. ZEUS ||140|| set a preliminary upper 



limit on the anomalous magnetic coupling of Ktu-y < 0.22 at 95%C.L. 

HI | ]141|| performed a search for single top production in e~^p (82 pb^^) and e^p data 
(14 pb~^) looking at events with either an isolated lepton plus missing transverse 
momentum, or hadronic events with three jets of large transverse momentum. Un 
upper limit of Ktu-y < 0.25 at 95%C.L. was derived. 

The HERA limits restrict strongly the allowed space for anomalous couplings of the 
top as shown in Fig. |5^ where the limits on the anomalous magnetic and vector 



couplings tuZ: Ja^ + are shown from HERA, LEP and the TEVATRON. 



6.7 Events with isolated leptons and large missing transverse 
momentum 

Events produced in ep collisions with high energy leptons and large missing transverse 
momentum are promising candidates for new physics such as production of squarks 
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Figure 53: NC cross sections da/dQ"^ normalised to the SM expectation for e~^p and 
e~p data. The curves show the effects of graviton exchange in large extra dimensions 
given by a common fit to the scale Ms with couplings A = +1 and A = — 1. The errors 
represent statistics and uncorrelated experimental systematics. The values shown for 
Ms represent 95% CL. lower limits. 
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Figure 54: Limits on the anomalous tM7 (magnetic) coupling k^^^ and on the anoma- 
lous tuZ (vector) coupling \ja\^v\ obtained at HERA, LEP and the TEVATRON. 



resulting from R-parity violating interactions. In the SM, events of this topology are 
predicted from W production where the W decays leptonically. However, since the 
cross section for W production is only about 1 pb, this contribution may constitute 
a small and manageable background. 

From a study of e^p interactions with 37 pb^^ HI has reported six events with an 



isolated, energetic lepton plus large missing transverse momentum [|142|] . For one of 
these events the lepton is an electron, for the remaining five it is a muon. The number 
of muons is significantly above the SM expectations. A similar search performed by 
ZEUS found agreement with the SM ||143|| . 



Preliminary results obtained from substantially larger data sets (HI: e+p with 82 
pb~^, e~p with 14 pb~^; ZEUS: e"*"p with 66 pb~^, e~p with 16 pb~^) have been 
presented in [p.44| , |145|| . 



For HI, the increased data set showed again a substantial excess over the SM predic- 
tion while ZEUS still observed good agreement with SM. For HI events the transverse 
momentum, p^, of the hadron system is plotted against the mass of the system 
formed by the lepton and the missing 4-momentum vector in Fig. |5^. The points 
with error bars show the observed events while the dots indicate the expected event 
distribution from W production which was calculated for a luminosity which is 500 
times larger than that for the data. In the muon sample several events stand out in 
a region where the contribution expected from W production is very small. 
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Figure 55: Distribution of events with an isolated lepton and large missing trans- 
verse momentum as a function of the transverse momentum of the hadronic sys- 
tem, and of the mass M^'^ of the system formed by the lepton and the missing 
4-momentum vector. The small dots represent the SM prediction for W production 
assuming a luminosity 500 times of the data sample. From HI. 
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In an effort to compare the HI and ZEUS data in a more transparent way, both exper- 
iments made an attempt to apply the same or at least similar selection criteria when 
selecting the topology of an isolated, energetic lepton plus large missing transverse 
momentum. In particular, HI restricted the track angular region to 0.3 < 9 < 2.0 
which is the range used by ZEUS. ZEUS required for the transverse momentum of 
the hadron system pf > 25 GeV. For the muon sample the missing transverse mo- 
mentum was required to be ptlmiss) > 12 GeV. Furthermore, events with a second 
/X of > 1 GeV were rejected. This requirement suppresses fi pair production by 
two-photon scattering, ep —>■ eNfifi, where is a proton or a low mass nucleonic 
state. None of the fi candidate events from ZEUS survived the latter cut. 

The number of observed and predicted events obtained with these criteria by the two 



experiments are presented in Table y |146|: 



Table 6: Observed and predicted (total and for VT-production alone) numbers for 
events with a lepton, transverse momentum pf of the hadron system above the cut 
values indicated and missing transverse momentum pt{'miss) > 12 GeV. From HI 
and ZEUS. 
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The numbers of events expected from SM processes as calculated by the two ex- 
periments agree well with each other. However, the numbers of events observed 
differ substantially: while HI sees a substantial excess over SM expectations, viz. 
for pf > 25 GeV 9 events are observed compared to 1.8 events predicted, the ZEUS 
event rates are in good agreement with SM: 1 event observed and 1.6 events expected. 

An increase in the data samples by factors 4-10 seems to be required before a 
definitive conclusion can be drawn. 



7 Outlook 



The year-2000 upgrade of HERA is expected to raise the luminosity by a factor 
of 3 - 5 which promises for HI and ZEUS each data samples as large as 1 fb~^. 
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At the same time the interaction regions for the two coUider experiments are being 
equiped with spin rotators in order to turn the transverse polarisation of electrons and 
positrons into a longitudinal polarisation. After the upgrade, the collider detectors 
will have improved vertexing capabilities for identifying c and b quarks with high 
efficiency. The combination of these improvements will increase substantially the 
physics capabilities ||147|| of the experiments. 



The charm and beauty part of the proton structure function can then be measured 
with high precision. A special feature of c and b quarks is their heavy mass which 
by itself provides a hard scale and which should suppress in deep- inelastic scattering 
possible contributions from nonperturbative processes. 

The structure functions F2, F^jxF^, F2, F2 can be measured between = 10 and 
40000 GeV^ in a single experiment. These measurements and the analyses of jet 
production in deep inelastic scattering are expected to determine the strong couphng 
with a precision of Aa^ = 0.0015, provided the necessary theoretical calculations are 
available in at least NNLO. 

With 1 fb~^ the study of diffractive processes in CC processes becomes a possibility. 

A luminosity of 1 fb~^ will allow to test quarks and electrons for substructure down 
to 4.10~^'^cm. The measured NC and CC cross sections will be sensitive to deviations 
from the Standard Model predictions down to the level of 30 - 50 MeV in M^. 
Additional W's and Z's can be detected up to 600 - 800 GeV. Perhaps most exciting 
will be the measurement of the charged current cross section for right (left) handed 
beam electrons (positrons). Any nonzero cross section contribution will be a sign for 
a new right handed charged current. 
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